Resonant excitation amidst dephasing: An exact analytic solution 
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An exact analytic solution is presented for coherent resonant excitation of a two-state quantum 
system driven by a time-dependent pulsed external field with a hyperbolic-secant shape in the 
presence of dephasing. Analytic results are derived for the amplitude and the phase shift of the 
damped Rabi oscillations. 
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I. INTRODUCTION 



Coherent resonant excitation represents an important 
notion in quantum mechanics [1, 2'|. Resonant pulses of 
specific pulse areas are widely used in a variety of fields in 
quantum physics, including nuclear magnetic resonance, 
coherent atomic excitation, quantum information, and 
others. Resonant excitation allows to establish a full con- 
trol over the quantum system, particularly in a two-state 
system, and realize any unitary transformation (qubit 
rotation) in it. Particularly important transformations 
include complete population transfer (transition proba- 
bility P = 1), complete population return {P = 0) and 
Hadamard transformation (P = ^). 

Crucial conditions for resonant coherent excitation are 
resonance (the frequency of the external field must be 
equal to the Bohr transition frequency) and coherence. 
Deviations from resonance (detuning) are detrimental 
and lead to rapid departure of the transition probabil- 
ity from the desired value. In this respect, an alternative 
to resonant excitation is provided by adiabatic excitation 
[3|, which is robust against such variations. 

Even more crucial for resonant excitation is coher- 
ence. Incoherent excitation, as described by Einstein's 
rate equations, allows only partial population transfer, 
e.g., at most 50% in a two- level system with equal de- 
generacies of the two levels fl| . Deviations from perfect 
coherence, which can be described by the quantum Liou- 
ville or Bloch equations, inevitably cause departure from 
the desired unitary transformation. Two general types 
of decoherence processes can be present: depopulation 
(e.g., due to spontaneous emission or ionisation) and de- 
phasing (e.g., due to elastic collisions, field fluctuations, 
coupling to the environment, etc.). 

In this paper, we present an exact analytic solution for 
pulsed resonant excitation of a two-state system in the 
presence of pure dephasing. Dephasing is recognised as 
one of the main obstacles in quantum information and the 
availability of precise analytic estimates of its effect can 
be very useful and important. We derive the exact solu- 
tion of the Bloch equation for a hyperbolic-secant pulse 
shape and a constant dephasing rate. We provide exam- 
ples of the general solution, which is expressed in terms 
of gamma functions, in various special cases of interest, 



e.g. for pulses with specific pulse areas. These results 
allow us to determine explicitly the deviations from the 
desired probabilities caused by the dephasing. 

This paper is organised as follows. In Sec. [Illwe de- 
scribe the model and present the exact analytic solution, 
which is derived in the appendix. In Sec. IIIII we analyse 
in some detail various special cases. In Sec. IIVI we give 
a summary of the results. 



II. ANALYTIC MODEL 
A. The Bloch equation 

Dephasing processes can be incorporated into the 
quantum-mechanical description of resonant excitation 
by including a phenomenological dephasing rate F into 
the Bloch equation fl|, y] , 



d 
dt 



u{t) 
v{t) 
w{t) 



-r -A(t) 

A(t) -r -n{t) 

n{t) 



u{t) 
v{t) 
wit) 



(1) 



This dephasing rate is the inverse of the transverse re- 
laxation time T2, r = I/T2 [1, 2]. Here u{t) = 2^pi2{t) 
and v{t) = 2Qpi2(t) are the coherences and w{t) = 
P22{t) — Pii(i) is the population inversion between the 
two states |1) and |2), with pmn {m,n = 1,2) being the 
density matrix elements. To be specific, we shall use the 
language of laser-atom interactions, although the results 
apply to any two-state system. In atomic excitation, the 
Rabi frequency r2(i) quantifies the time-dependent dipole 
interaction between the two states and it is proportional 
to the temporal envelope of the laser electric field E(t) 
and the transition dipole moment d, hn{t) — — d • E(t). 
The detuning A is the offset between the transition fre- 
quency of the two-state system loq and the laser field 
carrier frequency uj, A = loq — uj. 

We shall solve the Bloch equation ([1]) with the initial 
conditions 



u{—oo) ~ u(— 00) = 0, i(;(— 00) = — 1, 



(2) 



which correspond to a system initially in state |1): 
pii(— 00) — 1, p22(— 00) ~ 0. Our objective is to find 



the Bloch vector [u,v,w]'^ as i — > +00, and particularly, 
the population inversion w{oo), since the coherences van- 
ish at infinity due to the dephasing. 



B. The model 

We suppose an exact resonance, a hyperbolic-secant 
pulse and a constant dephasing rate, 



A{t) = 0, 

n{t) = f7oscch(i/r), 

r{t) = const. 



(3a) 
(3b) 
(3c) 



The dephasing rate F is a positive constant and T is the 
characteristic pulse width. The peak Rabi frequency Qq 
will be assumed also positive without loss of generality. 
The pulse area of the sech pulse (|3b|) is 



A 



+ OC 



n{t)dt = vrr^oT. 



(4) 



The exact solution in the coherent limit 



For F == 0, the Bloch equation ([T]) is solved exactly 



w{co) 



cos^. 



(5) 



Of particular interest are the cases when A is equal to 
an integer or half-integer multiple of tt. There are three 
cases of special significance. 

• The odd-i: pulses with area 

A^i2n + l)n, (n = 0,1,2,...), (6) 

invert the population, w{oo) = 1 (pn = 0, P22 = !)• 
A special case is the tt pulse with A — n. 



• The even-TT pulses with area 

A^2mr, (n = 0,1,2,...), 



(7) 



restore the population to the initial state, w{oo) — 
— 1 (pii — 1,P22 ~ 0). A special case is the 27r- 
pulse with A = 27r. 



• The half-integer-TT pulses with area 



A= {2n + l) 



(n = 0,1,2,...), 



(8) 



create an equal superposition between states 1 and 
2, zi;(oo) = {pii = P22 = 5)- A special case is the 
half-TT pulse with A — 7r/2. 

All these three cases are of great importance and such 
pulses are widely used in various applications in quantum 
physics, e.g. in nuclear magnetic resonance, coherent 
atomic excitation and quantum information. We shall 
therefore pay special attention to these cases in the ana- 
lytic solution, which we shall derive below. 



D. The exact solution with dephasing 



Because of the resonance condition (|3a| , it follows from 
Eq. Il]) that the equation for ii decouples (with the over- 
dot denoting a time derivative). 



u{t) = -Tu{t), 
and can be solved independently, 



u{t) = u{—oo)e 



-rt 



0, 



(9) 



(10) 



where we have used the initial conditions p|. 

We change the Bloch variable v{t) — ~ix{t) and Eq. 
(fT|) is reduced to two coupled equations, 

ix{t) ^ -iVxit) + n{t)w{t), (11a) 

iw{t) = nit)x{t). (lib) 

These equations resemble the Schrodinger equation for 
the Rosen-Zener model j^] with irreversible loss from one 
of the states [a, Q. For reader's convenience, the deriva- 
tion is adopted to our case and given in the Appendix. 
The exact solution for w{t) as t ^ -(-00 reads 



w(oo) 



rMI + 7) 



F(i+7 + a)F(i+7-a) 



(12) 



where T{z) is the gamma function [3, M, @l ^"^"^ the di- 
mensionless parameters a and 7 are defined as 



a = ^IqT, 



TT 



7 



(13) 



Because of the dephasing, the coherences vanish, 

u(oo) = v{(X)) = 0. (14) 

By using the reflection formula F(z)F(l — z) — 
7r/sin(7rz) 8], Eq. (fT2|) can be written also as 

wioo) = - ^' (^ + j| r (1 - 7 + <^) _ 

^ ' 7rF(i+7-fa) ^ /7 V y 



For 7 = 0, Eq. (|T5|) reduces to the lossless solution 
([5|). Equation (fT5|) shows that w(oo) vanishes whenever 
the cosine factor vanishes, i.e. for |a — 7I = ^7 1, f , ■ ' • ■ 
Hence the values of the pulse area, for which an equal su- 
perposition between states |1) and |2) is created (pn — 
P22 = 5)7 are shifted from their half-7r values ([8|), 



A={2n+1 + TT)-. 



(16) 



Several approximations to w{oo) are given below in 
some special cases. 
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FIG. 1: The population inversion, Eq. (|17|l for integer-vr 
pulses and Eq. (|20p for half-integer- tt pulses, against the de- 
phasing rate for different values of the pulse area A = ttq = 
7rf2oT, denoted on the respective curve. 

III. SPECIAL CASES 

A. Specific pulse areas 

For a — n (ni: pulse), where n is an integer, Eq. (|12p 
reduces to 



w{oo) 



n 

fe=0 



27 - 1 - 2/c 

27 + l + 2fc' 



(17) 



where we have used repeatedly the recurrence relation 
r{z + l) = zT{z). It follows from Eq. ^ that w(oo) = 
for 7 = ^ for any integer a = n ^ 1. 



For a = 1 (tt pulse) we have 

1-27 



w{oo) = 



1 + 27' 



(18) 



Hence the population inversion is a decreasing function 
of 7, which decreases to w^ = 1 — e for 



11 



7e 



We 



21 



(19) 



For ujj = 0.9, 0.5 and we find 7^ = g^i ^ and i, i.e. 
the inversion decreases very rapidly as 7 increases. 

Similar simple expressions can be derived from Eq. 
p?| for other cases of integer- tt pulses. 



For a - 
dni) that 



i (half-integer- TT pulse) we find from Eq. 



"^°°^ = ' rMi + 7) M^^' 



(20) 



The factor in front of the product gives u;(oo) for a = ^ 
(7r/2 pulse). 

In Fig. [1] the population inversion w{oo) is plotted 
against the dephasing rate for different pulse areas. 



B. Weak dephasing 

When 7 ^ 1 we use the relation r(a + 7) = 
r(a) [1 -|- '0(a)7 -|- 0(7^)] , where ^{a) is the psi func- 
tion, and the relation ip{l/2) = — (c-|-21n2), where 
c = 0.5772156649 ... is the Euler-Mascheroni constant 
0,[|. We thus find from Eq. ^ that 



w{oo) 



1-2 

X cos [tt {a — 7)] 



c + 21n2 + ^( - + a 



7 + 0(7') 

(21) 



The first factor {•••} describes the amplitude of the 
damped Rabi oscillations and the cos factor describes 
the phase of the oscillations. The maxima and the min- 
ima of these oscillations are shifted by 7r7 (if the small 
additional shift from the damped amplitude is neglected) 
from their coherent values © and ([7]), respectively. The 
factor {• • • } in Eq. ((2T|) displays explicitly the damp- 
ing of the amplitude and its departure from 1 as 7 rises 
from zero. Since ■0 (5 + a) is an increasing function of a 
[8|, the damping eifect is stronger for larger pulse areas, 
which is shown explicitly below. 

For a — J + n, near the nth extremum, we find from 
Eq. dn]) by using Eq. (6.3.4) of .8^ for V(n + 5) that 



w(oo) « (-1)"^ 



n 

' ^ 2A: - 1 



-0(7') 



fe=i 



n 


= 1 


-4, 




n 


= 1 




n 


= 2 




n 


= 3 




n 


= 4 



22D gives 



(22) 



(23a) 



wioo) w 1-47 + 0(7^), 

16 

w{oo)^-l + —-f + Oi-f^), (23b) 

92 

^(00)^1 7 + 0(7^), (23c) 

15 

u;(oo)« -1 + ^7-^0(7'). (23d) 



The cases of n = 1 and 3 correspond to the first and sec- 
ond maxima (tt and 37r pulses), and n = 2 and 4 to the 
first and second minima (27r and 47r pulses). Equations 
(|22| and (|23|) show explicitly how the values of the pop- 
ulation inversion for these mr pulses depart from their 
values ±1 as 7 departs from zero. These equations also 
demonstrate that the effect of dephasing is stronger for 
larger pulse areas (since the coefficient in front of 7 in- 
creases with n), which is indeed seen in Fig. [1] 

In Fig. [2] the population inversion is plotted against 
the pulse area for several values of the dephasing rate. 
As predicted, for 7 == ^ (FT = 1), the nodes of w are 
situated at pulse areas A = tt-tt, where in the absence of 
dephasing one finds the extrema. For 7 = 1 (FT — 2), 
the maxima are situated approximately at a = 2, 4, 6, . . ., 
where one finds the minima (even-7r pulses) for 7 = 0; 
likewise, the minima are situated approximately at a = 
3, 5,7,.. ., where one finds the maxima (odd-7r pulses) for 
7 = 0. 
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FIG. 2: The population inversion p2|) against the pulse area 
A = na for difterent values of the dephasing rate, FT = 
0.01, 0.1, 0.2, 0.5, 1, 2 (denoted on the respective curve). 

C. Strong dephasing 

For 7 ^ a, 1, the population inversion w{oo) has the 
asymptotics 

^,2 



w{oo) 



■ exp 



-- + 0(7-^) 

7 



(24) 



which is obtained from Eq. (|12p by using the Stirling 
asymptotic expansion [8], 



r(z) 



2-Kz'- 



' exp 



-2; + 



ifc + ^(Ni 



,(25) 



]argz| <7r, |z|>l). 



The inversion w{po) decreases in a Gaussian fashion 
against a. For large 7, w tends to its initial value of 
— 1, rather than to the incoherent limit it; = 0, which is 
a result of quantum overdamping 'G]. This behavior is 
indeed seen in Fig. [1] where we have verified that Eq. 
([24ll describes very accurately the asymptotic decrease of 
w(cxd) for large 7 (not shown for simplicity). 



D. Large pulse area 

When a 3> 7,1, the population inversion w(co) has 
the following behavior 



w(po) 



--r^ (^+7) oT'^^ [1 + 0(^-2)] cos7r(a-7), 

(26) 
which is derived from Eq. (jlSp by using Eq. (j25p . Equa- 
tion (p6)) shows that as a increases, the oscillation am- 
plitude vanishes as a"^'' and for sufficiently large pulse 
areas. 



A>,A, 



-KE 



rM^ + 7) 



(27) 



the population inversion decreases below e (|w| < e), i.e. 
the two-state system evolves towards a completely inco- 
herent superposition of states |1) and |2) (pn — p2i — \, 
P12 — 0). For instance, ^0.1 ~ 1-58 x IO^tt, 4157r, and 
3.187r for 7 = 0.1, 0.3, and 1, respectively. 



IV. CONCLUSIONS 

In this paper we have presented an exact analytic so- 
lution for resonant excitation induced by a pulse with 
a hyperbolic-secant shape in the presence of dephasing 
processes. The exact solution p2|) is given in terms of 
gamma functions. Dephasing affects the Rabi oscilla- 
tions in two ways: shifting the oscillation phase by ap- 
proximately -kTT /2 and damping the oscillation ampli- 
tude: the larger the pulse area, the stronger the damping. 
The implication is that one cannot reduce the dephasing- 
induced losses of efficiency by increasing the intensity of 
the field (e.g. replacing a tt pulse by a Stt pulse) since 
this will actually increase the losses. 

Various special cases of pulses with specific areas have 
been considered and various limits have been derived in 
terms of elementary functions. The results provide ex- 
plicit and simple estimates of the effect of dephasing on 
resonant excitation, e.g. in the cases of tt, 27r and 7r/2 
pulses, which are of great importance and widely used in 
many fields. 
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APPENDIX A: EXACT SOLUTION 

The first step in solving Eqs. (fTTjl is to decouple them 
by differentiating the equation for w and replacing v and 
-i), found from Eqs. (fTTjl : this gives 



*-r + ^ 



fi^w^O, 



(Al) 



with an overdot denoting d/ dt. We change the indepen- 
dent variable from t to z{t) — i[tanh(t/T) -f 1]; hence 
z(— 00) = and z{+oo) — 1. Then 

z{l-z)W" +(]-+"i-z\w' + a^W = Q, (A2) 

where ' = d/ dz, W [z{t)\ = 'w{t) and a and 7 are defined 
by Eqs. P^ . This equation has the same form as the 
Gauss hypergeometric equation [8, 9] , 

z(l - z)F" + [v~{\ + p + l)z] F' - XfiF = 0, (A3) 



upon the identification 
X = a, ^ 



(A4) 



The complete sohition of this equation, expressed by 
a superposition of two linearly independent solutions of 
Eq. (|A2|) , depends upon the value of ly. 

The case v ^ 1, 2, 3, .. .According to Sec. 9.153.1 of 
Q, the solution of Eq. (|A2p can be expressed in terms 
of the Gauss hypergeometric function [^, Q as 

W{z) = A^F{\ii-v-z) 

+A2Z^-''F{X + 1 - t/, ^ + 1 - z/; 2 - I/; ^)V5) 

From here and using Eq. (|lip . it can be found that 

V{z) 



iV^ 



A/i 



Ai-^F(A + 1, /I + 1; // + 1; z) 



+ ^2(1 - v)z-''F{\ + l-v,^i+l-v;l- vi^i) 



with V[z{t)] = v{t), where Eqs. (15.2.1) and (15.2.4) of 
[8| have been used. The integration constants Ai and A2 
can be determined from the initial conditions ((2|). 



Ai 



A2-O. 



Hence w{oo) — W(l) = -F(A, fi; v\ 1) or 



w{oo) 



T[u)T{y-X-y:) 
' T(y - X)T{v - ^lY 



(A7) 



(A8) 



where Eq. (15.1.20) of [Sj] has been used. Referring to 
Eqs. (HH), one obtains Eq. P^ . 

Equation (jASp has been derived under the assumption 
that z^ 7^ 1, 2, 3, . . .; then the two terms in Eq. (jA5p are 
linearly independent. Suppose now that v = n where 
n = 1,2,3,..., that is 7 = i,|,|,...; Then the two 



terms in Eq. (|A5|) are linearly dependent for t/ = 1 while 
the second term is not defined for v = 2,3,4, .. . 

The case i/ = 1. According to Sec. 9.153.2 of [^, the 
solution of Eq. (|A2p for i^ = 1 is 



W{z) = AiF{X,fi;l;z) 
+ A2 



FiX,,;l;z)lnz + ±^-^^z\^9) 



fc=i 



{kiy 



with (x)fe = r(a; + fc)/r(a;) and -01 = VXA + fc) - '0(A) + 
■0(m + k)- ^{fj.) - 20(fc + 1) + 2V'(1), -0(2;) being the psi- 
function [8|. Since the second term diverges for z = 0, 
the initial conditions ([2]) require Eqs. (|A7p to be satisfied 
and Eq. (|A8p applies again. 

r/ie case v = n + 1 (n — 1,2,3,...) and X,[i ^ 
0, 1, 2, . . . , n - 1. According to Sec. 9.153.3 of |9| and 
Eq. (15.5.19) of [i|, the solution of Eq. ((X2|) in this case 
is 



W[z) = AiF(A,/^;n + l;z) + A2 



F(X, ^i;n + I; z) Inz 



E 



{X)k{fJ.)k'>Pn^ 



{n + l)kk\ 



iz^ 



(fc-l)!(-n)fc 



{Alfe) 



,.^1 v^ . -,.-■- ^^(1-A)fe(l-M)fc 

with ■0„+i = '0'(A + fc) - ?A(A) + ip{^ + fc) - '0'(^) - ■0(ri + 
1 + fc) + V'(?^ + 1) - ^(1 + k) + ?A(1). Again, the second 
term diverges for z — and the initial conditions ^ 
require Eqs. (|A7p to be satisfied and hence, Eq. (jASp 
holds again. 

The case v = n + 1 (n = 1,2,3,...) anrf X or fi = 
0, 1, 2, . . . , n — 1. Suppose that A = m + l<n + l. Then, 
according to Sec. 9.153.4 of Q, the solution is given by 
Eq. (|A5p . in which the second hypergeometric function 
reduces to a polynomial in z^^. Since it diverges for 
z = 0, the initial conditions ^ require Eqs. (|A7p to be 
satisfied and Eq. (jASp holds again. 

In conclusion, for the model ([3]), in all cases the final 
population inversion w{+oo) is given by Eq. (jASp . 
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